UNCLASSIFIED 


Defense  Technical  Information  Center 
Compilation  Part  Notice 

ADPO 11621 

TITLE:  Magnetization  and  Giant  Magnetoresistance  of  the  System  of 
Interacting  Fine  Particles 

DISTRIBUTION:  Approved  for  public  release,  distribution  unlimited 


This  paper  is  part  of  the  following  report: 

TITLE:  International  Conference  on  Electromagnetics  of  Complex  Media 
[8th],  Held  in  Lisbon,  Portugal  on  27-29  September  2000.  Bianisotropics 
2000 

To  order  the  complete  compilation  report,  use:  ADA398724 

The  component  part  is  provided  here  to  allow  users  access  to  individually  authored  sections 
of  proceedings,  annals,  symposia,  etc.  However,  the  component  should  be  considered  within 
the  context  of  the  overall  compilation  report  and  not  as  a stand-alone  technical  report. 

The  following  component  part  numbers  comprise  the  compilation  report: 

ADP01 1588  thru  ADP01 1680 


UNCLASSIFIED 


143 


Magnetization  and  Giant  Magnetoresistance 
of  the  System  of  Interacting  Fine  Particles 

C.  Xu1,  Z.  Y.  Li2’1  and  I.  E.  Dikshtein3 


1 Department  of  Physics,  Suzhou  University,  Suzhou,  215006,  China 

2 CCAST(World  Laboratory),  P.O.Box  8730,  Beijing  100080,  China 

3 Institute  of  Radioengineering  and  Electronics,  Russian  Academy  of  Sciences, 

103907  Moscow,  Russian  Federation 

Fax:  +7-095-2038414;  email:  diksh@mail.cplire.ru 

Abstract 

The  magnetization  and  giant  magnetoresistance  (GMR)  of  nanosized  magnetic  particles  em- 
bedded in  a nonmagnetic  metallic  matrix  are  numerically  investigated.  By  considering  the 
classical  dipolar  interactions,  we  apply  a Monte-Carlo  simulation  technique  to  calculate  M 
vs  H,  GMR  vs  H,  and  GMR  vs  ( M/M„ ) with  M,  the  average  magnetization  along  the  nor- 
mal to  the  layer,  Ms,  the  saturated  magnetization,  and  H,  the  applied  magnetic  field.  It  is 
shown  that  the  interfacial  spin-dependent  scattering  of  conductance  electrons  is  dominant  in 
GMR  effect  and  the  distance  between  the  neighboring  particles  is  an  important  parameter 
to  obtain  the  GMR  effect,  while  the  size  distribution  only  modify  the  shape  of  the  curve  of 
GMR  versus  H. 

1.  Introduction 

The  discovery  of  the  giant  magnetoresistance  (GMR)  effect  in  inhomogeneous  alloys  of  magnetic 
and  nonmagnetic  metals  [1]  has  attracted  a great  deal  of  attention  to  these  materials.  They 
consist  of  nanosized  particles  or  clusters  (e.g.,  Co,  Fe,  Ni)  embedded  in  a nonmagnetic  metallic 
matrix  (typically  Cu,  Ag).  The  magnetic  transport  properties  of  granular  metals  are  concerned 
with  the  size  and  spatial  distributions  of  the  fine  particles  or  clusters  and  the  interaction  between 
the  particles.  Previous  works  studied  the  dependence  of  GMR  on  the  size  distributions  [2]  and 
successfully  explained  some  experimental  results.  The  interactions  between  particles  can  have 
a dipolar,  Ruderman-Kittel-Kasuya-Yosida  (RKKY),  or  a superexchange  character,  depending 
on  the  magnetic  properties  of  the  matrix.  Altbir  et  al.  [3]  found  that  the  classical  dipolar 
interactions  are  dominant  in  Co-Cu  systems. 

2.  Theory 

In  the  present  work  we  study  the  magnetization  and  GMR  effect  of  the  assembly  of  single- 
domain spherical  ferromagnetic  particles,  Each  particle  is  a saturated  single  domain  and  its 
magnetic  moment  pi,  and  the  direction  of  its  uniaxial  anisotropy  axis  is  random  in  space.  The 
particles  are  placed  in  a square  array  consisting  of  12x12  cells.  The  diameter  of  the  particle  i 
is  di  and  the  distance  between  two  neighboring  particles  is  ro-  After  considering  the  classical 
dipolar  interaction  and  crystalline  anisotropy  energy,  the  total  energy  of  the  system  for  a given 
configuration  {pi}  of  the  magnetic  moments  is 
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E{im)  = SE  ^ sin2  «<  - Pi  ■ H],  (1) 

i=l  j>i 

with  Eij , the  energy  of  the  dipolar  interaction,  Ku,  the  effective  anisotropy  constant,  Vi,  the 
volume  of  a particle  i , and  ctj,  the  angle  between  the  direction  of  the  crystalline  axis  and  pi- 
For  a given  temperature,  the  reduced  magnetization  < m > can  be  calculated  by  averaging 
m = M/Ms  = | Xu  fii\/ {Np)  over  cluster  configurations  after  thermal  equilibrium  has  been 
reached.  The  crucial  factor  for  GMR  in  granular  system  is  the  average  value  < cos  6ij  >,  where 
Oij  is  the  relative  angle  between  the  magnetic  moments  in  sites  i and  j [4],  It  implies  that  the 
magnetic  transport  properties  are  primarily  caused  by  the  spin-dependent  scattering  process 
of  conductance  electrons  between  magnetic  particles.  Namely,  the  spin-dependent  scattering  at 
interfaces  between  the  magnetic  particles  gives  rise  to  the  GMR  to  a greater  extent  than  the 
scattering  within  the  magnetic  particles  would  do  [5].  In  the  case  that  the  distance  between  the 
neighboring  particles  does  not  far  exceed  the  electronic  mean  free  path  A [6],  the  variation  in 
resistivity  of  a granular  system  with  the  degree  of  field-induced  magnetic  order  may  be  simply 
pictured  as  p = po  — k < cos  >,  where  po  and  k are  constants.  Assuming  that  there  are 
no  correlations  between  the  magnetic  moments  of  particles,  the  magnetoresistance  A p/p  can  be 
written  as  A p/p  — —(k/po)  < cos  o\^  >2=  — ( k/po)m 2.  Such  a quadratic  dependence  of  A p/p 
on  m is  actually  found  by  some  experiments  [1].  However,  other  experiments  [2]  showed  that 
A p/p  does  not  vary  quadratically  with  m because  of  the  size  distribution  of  magnetic  particles 
and  interaction  between  them. 

For  the  system  involving  coupling  between  magnetic  particles,  the  assumption,  that  the  av- 
erage value  < cos  >o  for  H = 0 is  not  vanished,  gives 

Ap/p  = (<  cos  e\f  >0  - < cos  o\?  >)(Q~  < cos  >0),  (2) 

where  Q — po/k  is  the  field-independent  constant. 

The  thermal  averages  of  the  system  above  are  obtained  using  the  standard  MC  procedure  and 
the  Metropolis  algorithm  [7].  The  system  is  assumed  to  have  reached  thermal  equilibrium  after 
104 *  Monte-Carlo  steps  per  spin.  Then,  we  are  able  to  get  the  thermal  averages  as  an  arithmetic 
average  over  the  accepted  configurations  (500  accepted  configurations  for  ensemble  averages), 
and  to  calculate  the  < m > and  GMR.  Data  for  our  MC  simulation  is  generated  and  calculated 
as  follows.  Each  particle  is  assigned  a random  crystalline  anisotropy  ( Ku  = 4.0  x 106 erg/cm3) 
and  a random  direction  of  magnetic  moment  at  initial  state.  These  particles  are  placed  in  the 
magnetic  field  H applied  along  the  normal  to  the  array.  The  distance  ro  between  the  particles 
was  taken  as  6.0  nm  (except  for  the  up  triangles  in  Fig.  1)  which  is  comparable  in  magnitude  with 
A [6].  Single-domain  ferromagnetic  particles  exhibit  the  phenomenon  of  superparamagnetism 
and  the  blocking  temperature  for  H = 0 of  a particle  of  diameter  d— 4 nm  is  equal  to  38  K 
[8].  So  we  choose  T=40  K close  to  Tb(0). 

3.  Results 

First,  for  simplicity,  we  choose  the  same  diameter  dj(=3  nm)  of  all  the  magnetic  particles,  which 
is  the  typical  average  size  of  particles  for  granular  materials  [6].  In  Fig.  1 we  plot  the  graph  of 
magnetization  M/Mq  vs  H.  Four  different  sets  of  data  are  shown  that  correspond  to  a system 
of  particles  with  random  anisotropy  only  (squares;  ro  = 6 nm),  a moderate  dipolar  system 
with  random  anisotropy  (circles;  ro  = 6 nm),  a system  with  moderate  dipolar  interaction  only 

(triangles  down;  Ku  = 0,  ro  = 6 nm),  and  a strong  dipolar  system  with  random  anisotropy 

(triangles  up;  ro  = 3 nm).  We  notice  that  M exhibits  the  different  field  dependence,  depending 
on  the  interplay  of  the  single-particle  anisotropy  and  the  dipolar  interaction  effects.  E.g.,  for  a 
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Fig.1 . The  dependence  of  M/Mt  on  H for  T=40  K.  The  symbles  are 
obtained  when  we  do  not  consider  the  dipolar  interaction(DI) 
(squares),  or  do  not  include  the  crystalline  anisotropy  (CA) 
(down  triangles).  The  circles  and  up  triangles  are  obtained  by 
considering  both  the  D1  and  CA  when  the  r equals  6nm  and 
3nm  respectively. 


Fig.2.  The  dependenc  of  GMR  on  H for  T=40  K.  The  symbles  are  obtained  when 
r^=6nm,  d=3nm  (circles),  d=2-4nm(squares),  d=0~6nm{  down  triangles),  and 
r=3nm,  d=3nm  (up  triangles). 


system  with  random  anisotropy  only  the  magnetization  reversal  of  the  assembly  of  particles  is 
well  described  by  both  a coherent  rotation  of  pi  from  the  easy  axes  to  the  direction  of  H and 
by  a thermally  activated  process  over  the  anisotropy  barrier. 

A pure  dipolar  system  ( Ku  = 0)  For  H = 0 possesses  the  interaction-induced  anisotropy. 
Out-of-plane  orientation  M has  a large  energy  due  to  the  demagnetization  effect.  Therefore 
in-plane  configuration  of  pi  is  realized.  For  the  free  boundary  conditions  accepted  here  the 
ground  state  is  antifferromagnetic  due  to  the  demagnetization  effect  of  the  lateral  boundaries. 
This  effect  is  not  so  strong  compared  with  the  demagnetization  effect  of  the  surface.  Once  again, 
if  a field  is  applied  to  the  array  the  magnetic  moments  pi  rotate  coherent  to  the  direction  of  H. 

If  both  the  dipolar  interaction  and  the  random  uniaxial  anisotropy  of  particles  are  involved  it 
is  more  difficult  for  M to  reach  saturation  since  these  effects  impede  the  collinear  ordering  of  pi 
and  reduce  the  magnetization  of  the  system.  The  anisotropy  induced  by  dipolar  interactions  is 
very  sensitive  to  the  spatial  arrangement  of  the  particles  (average  distance  between  the  particles 
and  their  size)  and  it  is  enhanced  with  decreasing  ro . Therefore  in  a strong  dipolar  system 
(triangles  up;  ro  = 3 nm)  the  M is  harder  to  be  saturated  in  comparison  with  a moderate 
dipolar  system  (circles;  ro  = 6 nm). 

In  Fig.  2 we  plot  the  field  dependence  of  A p/p  for  T=40  K.  All  curves  are  calculated  in  terms 
of  Eq.  (2)  with  Q=5.2  after  considering  both  the  effects  of  anisotropy  and  interactions  when 
the  system  reaches  the  equilibrium  state.  To  investigate  the  influence  of  the  distance  ro  and  the 
particle-size  distribution  on  the  GMR  four  different  sets  of  data  are  shown  that  correspond  to 
particles  of  fixed  diameter  with  strong  (up  triangles;  dj=3  nm,  ro=3  nm)  and  weak  (circles;  d;=3 
nm,  ro=6  nm)  dipolar  interactions,  a narrow  particle-size  distribution  from  2 to  4 nm  (squares; 
d—  3 nm,  ro=6  nm),  and  a wide  particle-size  distribution  from  0 to  6 nm  (down  triangles;  d=  3 
nm,  ro=6  nm)  with  d is  the  average  particle  size.  From  Fig.  2 it  follows  that  the  GMR  effect 
depends  crucially  on  the  particle  density.  If  the  dipolar  interaction  are  strong  (</,— 3 nm,  ro=3 
nm)  all  the  particles  are  coupled  to  each  other,  and  the  negative  GMR  disappears.  For  a narrow 
particle-size  distribution  from  2 to  4 nm  (squares)  the  field  dependence  of  A p/p  is  closed  to 
that  for  the  system  of  particles  with  fixed  diameter  and  a weak  dipolar  interection  (circles).  In 
this  case  the  blocking  temperature  of  the  maximum  particle  Tj,  «38  K,  and  all  the  particles  are 
superparamagnetic.  This  result  was  observed  in  [5].  For  a wide  particle-size  distribution  from 
0 to  6 nm,  certain  of  large  particles  are  blocked  at  40  K in  the  region  of  strong  magnetic  fields, 
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whereas  the  rest  small  particles  are  still  in  the  superparamagnetic  state.  Then  the  modulus  of 


\Ap\/p  is  substantially  reduced  relative  to  the  case 
of  a narrow  particle-size  distribution.  This  effect 
was  observed  and  explained  in  [2]  for  the  systems 
with  a wide  particle-size  distribution.  So  we  be- 
lieve that  the  size  distribution  strongly  affects  the 
shape  of  the  curves  of  GMR  vs  H. 

In  order  to  investigate  in  more  detail  the  role  of  a 
particle-size  distribution  in  GMR  effects,  we  show 
in  Fig.  3 the  dependence  of  A p/ p vs  M/Ms  for 
a narrow  (hollow  squares)  and  wide  (up  triangles) 
particle-size  distributions.  From  Fig.  3 it  follows 
that  for  a narrow  particle-size  distribution  from  2 to 
4 nm,  the  behavior  of  GMR  vs  ( M/Ms ) is  close  to 
the  parabolic  line,  whereas  for  a wide  particle-size 
distribution  from  0 to  6 nm,  it  is  deviated  essen- 
tially from  the  parabolic  law.  It  indicates  that  the 
wide  distribution  of  magnetic  particles  may  explain 
the  noncompliance  with  the  parabolic  law  for  the 
GMR  as  a function  of  the  M [2]. 


M/M, 


Fig. 3.  The  dependence  of  GMR  on  M/M4  for  T=40  K.  The  hollow 
squares  and  up  triangles  are  obtained  by  varing  of  r.  from 
2 to  4 nm  and  0 to  6 nm  respectively.  The  solid  line  fs  the 
parabolic  line. 


4.  Conclusion 

In  conclusion,  we  have  presented  the  results  for  the  field  dependence  of  M and  GMR  of  a granular 
magnetic  film.  We  demonstrate  that  the  magnetic  properties  of  the  system  depends  essentially 
on  the  particle-size  distribution  and  the  average  distance  between  the  magnetic  particles.  To 
discuss  experimental  data  the  state  of  single-domain  magnetic  particles  is  usually  assumed  to 
be  blocked  or  collective  at  low  temperatures  and  to  be  superparamagnetic  at  high  temperatures. 
Within  our  MC  approach,  there  is  no  need  for  making  o priori  assumptions  about  the  particle 
state.  However,  it  is  likely  that  the  high  density  regime  favours  the  collective  state,  and  the 
low  density,  wide  particle-size  distribution  and  strong  magnetic  field  regime  favours  the  blocked 
state.  The  collective  behaviour  at  high  particle  density  reveals  itself  in  the  disappearance  of  the 
negative  GMR  effect.  A manifestation  of  the  blocked  state  effects  is  a substantial  decrease  in 
\Ap\/p  for  the  low  density  and  wide  particle-size  distribution  in  the  region  of  strong  magnetic 
field. 
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